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Abstract: Shelterbelts are important in defending against natural disas¬ 
ter and maintaining ecological balances in farmland. Understanding of 
the shelterbelt vegetation fraction is fundamental to regional research of 
shelterbelts using remote sensing. We used SPOT5 imagery with 
1 Ox 10m spatial resolution in combination with knowledge of the charac¬ 
teristics of shelterbelts to develop a method for retrieval of the vegetation 
fraction of shelterbelts by the pixel un-mixing model. We then used the 
method to retrieve values for shelterbelts in study area. By combining the 
parameters of photographic images with characteristics of shelterbelts, 
we developed a method for measuring the vegetation fraction of shelter- 
belts based on an advanced photographic method. We then measured the 
actual values to validate the retrieval result. The multiple correlation 
coefficients between the retrieved and measured values were 0.715. Our 
retrieval and measuring methods presented in this paper accurately re¬ 
flect field conditions. We suggest that this method is useful to describe 
shelterbelt structure using remote sensing. 

Keywords: shelterbelt; vegetation fraction retrieval; vegetation fraction 
measuring; remote sensing 

Introduction 

Shelterbelts, linear arrays of trees and shrubs planted to create a 
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range of benefits, are a major category of agroforestry practices 
(Buck et al. 1999). They protect crop fields by reducing wind 
erosion and improving crop water use, and they increase crop 
yields and economic returns (Kort 1988). However, the devel¬ 
opment of shelterbelt construction makes it more difficult to 
manage the shelterbelts. To efficiently acquire shelterbelt infor¬ 
mation on regional scales, remote sensing (RS) had been widely 
applied in shelterbelt research, including shelterbelt extraction 
(Wiseman et al. 2009) and landscape assessment of shelterbelts 
(Zhou et al. 1994; Shi et al. 2011). But little research has been 
undertaken to describe shelterbelt structures by use of RS. 

The vegetation fraction, the ratio of vegetation occupying a 
unit area, is an important parameter in developing climate and 
ecology models. Accurately acquiring the vegetation fraction of 
shelterbelts is necessary to study shelterbelt structure and growth 
by RS. 

There are many methods for retrieval of the vegetation frac¬ 
tion, the spectral un-mixing model being the most popular. In 
this research we used the sub-pixel model combined with knowl¬ 
edge of shelterbelt characteristics and site conditions to improve 
key parameters of the model. We then developed a method suit¬ 
able for extracting the vegetation fraction of shelterbelts. 

There are many methods for measuring site vegetation frac¬ 
tions, one of which that is widely used is advanced digital pho¬ 
tography (White et al. 2000). In the present study we used this 
method in combination with knowledge of shelterbelt character¬ 
istics and the camera parameters to develop a method for meas¬ 
uring the vegetation fraction of shelterbelts. We then used meas¬ 
ured values to validate the retrieval values. 

Materials and methods 

Data source 

Shelterbelts are linear objects. By field measurement, the canopy 
width of shelterbelts in northeast China is about 15 m. So the 
SPOT5 image with 10 m x 10 m spatial resolution was chosen as 
our data source. Shelterbelts are distributed in farmland. The 
spectral features are similar between shelterbelts and crops, and 
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this affects the extraction of the shelterbelt from farmland during 
the growing season. Therefore, the best date to acquire satellite 
images is in May, when the spectral differences between shelter- 
belt and crop are greater. The acquiring date of SPOT5 imagery 
in this study was April 29, 2007. 

Site description 

The coverage area of one SPOT5 image was selected as our 
study area, which was located in a wind damage zone of the 
Songnei Plain. This is a semi-humid zone mainly overlapping 
Dehui and Nong’an counties of Jilin Province (Fig. 1). This area 
is important for commercial grain production and typical shelter- 
belt construction in China. The shelterbelt construction of this 
area started in the 1950s. After sixty years of construction and 
management of shelterbelts, especially since the Shelter Forest 
System Program in the Three North Regions of China in 1978, 
this area has basically cultivated shelterbelt networks in all suit¬ 
able farmlands. The shelterbelts are important in ensuring grain 
production. 
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Fig. 1 Location of the study area 

Retrieval of the vegetation fraction of shelterbelt 
Dimidiate pixel model 

The dimidiate pixel model is common in the spectral un-mixing 
model, and it is a simple and practical remote sensing model (Li 
et al. 2004; Tian et al. 2004; Chen et al. 2001). The model as¬ 
sumption is that the surface of a pixel is covered with vegetation 
and non-vegetation. The Normalized Difference Vegetation In¬ 
dex (NDVI) is a type of quantitative value calculated from spec¬ 
tral information of surface objects received from remote sensors 
and reflects the condition of surface vegetation. 

According to the dimidiate pixel model, one pixel of NDVI 
value can be marked as V veg from green vegetation information 
and V so ii from bare soil information. Hence, the remote sensing 
model of the estimation of f c is defined as: 


fc=(N-N^,)l(N mg -N mi 0 (i) 

where, V soil and V veg represent the NDVI values of the bare soil 
or non-vegetation coverage pixel and a pure vegetation pixel 
respectively. N is the NDVI value in the equation. Theoretically, 
the V soil value of most bare soil surfaces approaches zero. How¬ 
ever, the range of V soil values is between -0.1 and 0.2 due to the 
impact of many factors (Carlson et al. 1997; Bradley 2002). V veg 
represents the maximum value of the pure vegetation pixel. Nev¬ 
ertheless, V veg also changes with time and space due to the effect 
of different vegetation types. Therefore, the determination of 
V soil and V veg values is critical to the effectiveness of this model. 

Determination of N so n and V veg 

There are many methods to obtain end-member values. Scarth et 
al. (2000) matched the geometric-optical model output to field 
values, and determined suitable end-member values for use in the 
spectral un-mixing model. This method had relatively accurate 
precision, but it required field measurement to ground-truth the 
image. Zeng et al. (2007) calculated the end-member values with 
the support of extensive in situ measurements by field spec¬ 
trometer and Hyperion image. The field measured end-member 
values were accurate, but because of atmospheric influence and 
differences in spatial resolution, the values on the satellite image 
can have more errors than field measurement. Gutman et al. 
(1998) made the minimum value in the image as N soih and the 
maximum value as V veg . However, the image values were easily 
disturbed by shadows or clouds, resulting in abnormal values in 
images. Toby et al. (1997) suggested that the ideal V veg is lower 
than 0.05 as a maximum value. But this was only an experimen¬ 
tal value that would be changed for different images. By analyz¬ 
ing the cumulative probability distribution of NDVI about vege¬ 
tation and soil, Li et al. (2004) chose the NDVI value with confi¬ 
dence of 5% as V veg . 

Based on review of the methods mentioned above, the method 
to determinate end-member values in this study was as follows: 
We used land use and soil type data, extracted NDVI values of 
woodland and cropland for each soil type, and chose the NDVI 
value with given confidence as V veg and V soil . 

First, we overlapped land use and soil type data. For V soil , the 
land use type ‘cropland’ was taken as a mask to extract NDVI 
values for each soil type. Then the NDVI values were sorted in 
descending order. The value with confidence of 0.5% was cho¬ 
sen as V soil for each soil type (Table 1). For V veg , if V veg in each 
soil type was calculated, the samples of woodland in the study 
area might be inadequate. Therefore, the NDVI values of crop¬ 
land for the entire study area were extracted and sorted in as¬ 
cending order. The value with confidence of 0.5% was chosen as 
V veg , and this was 0.31 in this study. 


Table 1. The 7V soil of each type of soil 


Soil type 

Black Soil 

Chernozem 

Meadow Soil 

Alluvial Soil 

N soil 

0.08 

0.08 

0.06 

0.07 
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Calculation of vegetation fraction of shelterbelt 
A SPOT5 image with 1 Ox 10m spatial resolution was chosen as 
our data source. After geometrical correction and radiometric 
correction, the NDVI ( N) was calculated by formula (2): 

N = (b\-b2)l(b\ + bl) (2) 

where, b\ is the apparent reflectance of the near-infrared band, b 2 
is the apparent reflectance of the red band. 

Formula (1) was used to calculate the vegetation fraction of 
shelterbelts. Where, N veg was 0.31, and A soi i can be found in Ta¬ 
ble 1. The calculated result is shown in Fig. 2. 



L 


Fig. 2 The retrieval result of vegetation fraction of the shelterbelt 

Result and analysis 

Validation method 

To select a sampling interval required understanding of the area 
captured on a photograph. The parameters of the camera used in 
this study are as follows: focal length 35 mm, diagonal direction 
62°, horizontal direction 53°, and vertical direction 37°. We as¬ 
sumed tree crowns to be oval-shaped. The parameters of the tree 
are shown in Fig. 3. We took the minimal visual angle (vertical 
direction) as standard, the actual length of photo in vertical direc¬ 
tion was calculated as: 2 m x 12.5 mx (37°/2) ~ 8.36 m. 

When the digital camera method was used to measure the veg¬ 
etation fraction, the margin of the photo was distorted. To guar¬ 
antee the precision of measurements, it is better to remove the 
margins of photographs. We defined a circle centered on the 
same spot as a photograph and with a diameter of 5 m as the 
valid area. 

Considering of the spatial resolution of a SPOT5 image and 
the valid area of a digital photograph, we set the sample interval 
at 5 m in this study. In Fig. 4, the square shaded area represents 
the pixel of a satellite image and the circle is the valid area of a 
digital photograph. The vegetation fraction of each photo within 


the valid area was computed by threshold, and finally, the mean 
value of 3x3 vegetation fractions in valid areas was taken as the 
vegetation of a pixel. 



Fig. 3 Model of taking photo of tree 



Fig. 4 Sampling interval diagram for measuring the belt vegetation 
fraction 

Results 

We chose an entire shelterbelt in our study area as an object, and 
measured the vegetation fractions on 9 June 2010 using the vali¬ 
dation method described above. The measured values and re¬ 
trieval values are shown in Fig. 5. Absolute error was relatively 
large and means absolute error was 0.23. Large error was caused 
mainly by the temporal difference between sample dates: the 
measured value was acquired in June, and the retrieval value was 
acquired in late April. Fig. 6 shows that the multiple correlation 
coefficient between the two values was 0.715, indicating high 
correlation. Additionally, the fluctuating trends yielded by both 
methods were consistent, indicating that the retrieval value could 
reflect the actual vegetation fraction of the shelterbelt. 
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Fig. 5 Comparison between measured vegetation fraction and re¬ 
trieved value 
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Fig. 6 Correlation between measured vegetation fraction and re¬ 
trieved value 


Conclusions 

We used the dimidiate pixel model to retrieve the vegetation 
fraction of a shelterbelt using SPOT5 image. The key parameters 
- end-member values were determined using land use and soil 
type data. To validate the accuracy of the retrieval result, we 
developed a validation method by computing the valid area of a 
digital photograph to select a sampling interval. We then used 
the method to measure an entire shelterbelt in our study area. The 
retrieval and measurement methods for determining the vegeta¬ 
tion fraction of shelterbelts developed in this paper were practi¬ 
cable. The retrieval result reflected the actual value of the vege¬ 
tation fraction, and this method can provide technology for de¬ 
scribing shelterbelt structures through the use of RS. 
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